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Abstract 

Spintronics is the generic term that describes magnetic systems coupled to an electric generator, 
taking into account the spin attached to the charge carriers. For this topical review of Spin 
Caloritronics, we focus our attention on the study of irreversible processes occuring in spintronic 
devices, that involve both the spins of the conduction electrons and the ferromagnetic degrees of 
freedom. The aim of this report is to clarify the nature of the different kinds of power dissipated 
in metallic ferromagnets contacted to an electric generator, and to exploit it in the framework of 
the theory of mesoscopic non-equilibrium thermodynamics. The expression of the internal power 
(i.e. the internal entropy production multiplied by the temperature) dissipated by a generic system 
connected to different reservoirs, allows the corresponding kinetic equations to be derived with the 
introduction of the relevant phenomenological kinetic coefficients. After derivation of the kinetic 
equations for the ferromagnetic degrees of freedom (i.e. the Landau-Lifshitz equation) and the 
derivation of the kinetic equations for the spin-accumulation effects (within a two channel model) , 
the kinetic equations describing spin-transfer are obtained. Both spin-dependent relaxation (usual 
spin-accumulation) and spin-precession in quasi-ballistic regime (transverse spin-accumulation) are 
taken into account. The generalization of the Landau-Lifshitz equation to spin- accumulation is then 
performed with the introduction of two potential energy terms, that are experimentally accessible. 

PACS numbers: 75.40. Gb,72.25.Hg,75.47.De 
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I. SUMMARY 



The approach of spintronics adopted here is that of Non-Equilibrium Thermodynamics 
1-5] applied at the Mesoscopic scale jf], 7] (MNET). The analysis is based on the expres- 
sion of the entropy production, i.e. on the expression of the power dissipated through the 
different relaxation mechanisms that characterize the system. The theory is adapted to the 
description of a plurality of out-of-equilibrium sub-systems exchanging energy, in which the 
role of environmental microscopic degrees of freedom are reduced to transport coefficients 
(damping, diffusion coefficients, conductivity, thermoelectric power, gyromagnetic ratio, and 
other Onsager coefficients) for the collective variables under consideration. The description 
holds at the mesoscopic scales, under the hypothesis of local equilibrium extended to internal 
degrees of freedom. The aim of this report is to propose an application of MNET to metallic 
spintronic devices that includes uniform ferromagnetic degrees of freedom explicitly. 

Before presenting the detailed derivation of the kinetic equations in the forthcoming 
sections II, III and IV, let us first summerize the general scheme of the report. In this intro- 
ductory section, the usual thermoelectric effect is first presented, and extended to the case 
of a bi-valuated internal variables: this is the two-channel model. An analogous approach 
is then performed in the space of the ferromagnetic degrees of freedom S. The coupling be- 
tween the E space and the internal degree of freedom of the electronic sub-system leads us to 
the spin-transfer kinetic equations. The detailed treatment of the ferromagnetic transport 
is given in section II, the detailed treatment of the two channel model of electric transport 
is performed in section III, and the coupling between both subsystems, i.e. the spin-transfer 
, is presented in the last section. 



A. Thermoelectric effects 



The internal power dissipated by a system is given by the internal entropy production 
(where S is the entropy of the system) multiplied by the temperature T. In the case 
of electric charges moving in one dimension z in a wire of section unity (i.e. the wire is 
contacted to two reservoirs of electric charge) and maintained at uniform temperature , the 
power dissipated inside the system is given by the Joule heating, i.e. the product of the 
electric current by the electric field 8: 
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dt e dz y ' 

where J e is the electric current, e the absolut value of the electric charge, and /i e is the 

electrochemical potential. The electric field is given by £ = —\^r- 

The application of the second law of thermodynamics ^ > leads us to define a first 

positive Onsager coefficient a (which is a function of the state variables) in order to built a 

positive quadratic form. The electric current writes: 

r = - a - 8 -f (2) 

e dz 

which is Ohm's law and a is the electric conductivity. 

On the other hand, the power dissipated inside a wire contacted to two heat reservoirs is 
the product of the heat flow J Q and the conjugated force (^) multiplied by T: 

T <£ = T J*.<>a) (3) 
dt dz\TJ y ' 

The Fourier's law is deduced from the second law of thermodynamics after introducing a 
positive Onsager coefficient k: 

JO = -«f (4) 

where k is the thermal conductivity. 

If both electric and heat reservoirs are contacted to the same wire, the internal power 
dissipated is: 

T^-J-^Tjo'a) (5) 



dt e dz dz \T y 

Now, beyond Ohm's law and Fourier's law, the two currents are coupled through the 
relevant Onsager thermoelectric cross- coefficients, 

r = + So— 

e dz dz (f>\ 

TO STadfx e , ~ ,dT y ' 

J = -jr~( K + TS a ) 7T 

e dz v ' dz 

where the Onsager cross-coefficients are expressed with the help if the Seebeck coefficient 

S, defined at zero electric current by the relation £ = SdT/dz and the conductivities a 
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and k. Note that according to the Onsager reciprocity relations, the Peltier coefficient IT, 
defined without electric field by the relation Jq = UJ e , verify to the relation II = TS. The 
existence of the cross-coefficients is justified by the fact that the charge carriers are also 
contributing to the transport of heat, or inversely, the transport of heat is contributing to 
the transport of electric charges. If a detailed microscopic theory is possible in the case of 
the Ohm's law and the Fourier's law, this is no longer the case in general for thermoelectric 
effects or other cross-effects. However, the knowlege of the detailed mechanisms of heat 
transport due to electric carriers is not necessary in order to derive Eq. (J6j). This justifies 



the interest of a thermokinetic phenomenological appoache applied to spin ca 



which the underlaying relaxation mechanisms are also not well known 
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B. Two-channel relaxation and spin-accumulation 

Let us assume that the ensemble of electric charges is composed of two different popula- 
tions. The difference is introduced through an internal degree of freedom, restricted here to 
a bi-valuated variable that takes the value a and 7. In the context of semiconductor physics, 
the two channel model was introduced in order to describes the transport of both electrons 
and holes 0, Q . In the context of the usual spin-accumulation effect due to spin-flip relax- 
ation, the two channels account for the spin up or spin down attached to the conduction 



electrons: a =f an d 7 =1 12Hl5|. However, from the point of view adopted here, this 
scheme should be generalized to a band structure in order to account for the s — d relax- 
ation mechanisms, that are responsible for the coupling between the spin of the conduction 
electron (mainl y s electron band) and the ferromagnetic degrees of freedom (related to d 



electron band) [l6(. In this case, electronic transport is also spin-dependent because the d 
band is full for majority spins (e.g. f) in usual 3d metallic ferromagnets [17 ]. 

However, without entering into the complexity of the spin-dependent relaxation mech- 
anisms, it is easy to generalize Ohm's law with adding the parameters a and 7 to the 
transport coefficients. A third kinetic equation should be introduced in order to take into 
account the power dissipated by the a — > 7 relaxation (spin-flip relaxation or s — d relax- 
ation). This relaxation is formally equivalent to a chemical reaction, driven by the chemical 



affinity A/x = \x a — /x 7 15]. Indeed, the power dissipated by the system reads then: 
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T 



dS[ 
dt 



-j: 



a ' e dz r e dz 



(7) 



where the flux of particles relaxing from one channel to the other (relaxation occurring 
in the space of the internal degrees of freedom) is given by The corresponding kinetic 
equations are deduced, after introducing a supplementary Onsager coefficient L. 



V 



cr a dfi e a 

e dz 
cr 7 <9/J* 



(8) 



e dz 
if = LA/i 

The set of equations Eqs (jSJ) is sufficient and necessary in order to describe, in the station- 
ary regime, spin-accumulation effects or any non-equilibrium contribution to the resistance 
due to A/i occurring at an interface. The corresponding effects (spin accumulation, giant 
magnetoresistance, etc) will be discussed in Section III. Onsager cross-coefficients may also 
be added at this stage of the analysis [18 ] . 

It is important to note the introduction of the parameter A/i in the irreversible processes 
described in Eq. (EJ): in spintronics, A/i is called "spin accumulation" and it was first 
introduced in this context by Van Kempen et al. [13]. This parameter plays the role of the 
pumping force that is responsible for the out-of-equilibrium relaxation occurring from one 
channel to the other. It will be also responsible for spin-transfer as described at the end of 
this report. 

Accordingly, it is convenient to rewrite Eq. |8| as a function of the variable A/i. Let 
us define the conductivity asymmetry by the parameter j3 such that = aa ~^ 1 and the 
mean conductivity 2<7o = a a + o 1 . On the other hand, the spin-polarized electric current 
is 5J e = J e a — and the spin- independent current is Jq = + J^. In this new system of 
equations, the Onsager matrix re-writes: 



/ 



CT (3oq 

f3a (T 
L 



\ ( -1 dA^ \ 

e dz 



/ 



e dz 

A/i 



(9) 



where /xq = /i« + A* 7 - 



5 



The generalization of the thermoelectric effect to the two channel case is straightforward, 
assuming that the electrons are thermalized (i.e. the temperature is the same for each 
channel) : 



•J ^ 



(T a dfi% dT 

f; r <J a <J a —— 

e oz oz 



J% = -S a TaJ£. - ( Ka + TSla a ) ^ 



r 



J? 



dz 



(10) 



e dz 



-<S 7 Tcx 7 — v „ 7 , - „ r 

where Si and i = {0,7}, are respectively the Seebeck and Fourier coefficient of each 
channel. However, in this report, we will not further investigate this set of equations. Some 



consequences in relation with experiments have been investigated and reported 
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set of equations (FlDjl . or other theoretical descriptions beyond the two channel model 29 



he 



investigated the caloritronic properties of spintronic systems. In contrast, the goal of this 
report is to describe the transport properties of the ferromagnetic degrees of freedom and 
the consequences of its interaction with spin-dependent electric sub-systems. 



C. Introduction of the ferromagnetic degrees of freedom 

Let us introduce the ferromagnetic degrees of freedom. This observable is defined in a 
physical space that is not the usual space considered above, but a space of magnetic 
moments S. This space should be considered as a space of internal degrees of freedom, in 
the same sense as the bi-valuated variable {a, 7 } of the two channel model (and in the same 
sens as the spin space in quantum mechanics) . The power dissipated by the magnetic system 
is then given by the flux Jq of magnetic moments (defined in the corresponding vectorial 
space E) multiplied by the magnetic force: 

T^ = J?.V E ^ (11) 

where Hq is the ferromagnetic chemical potential and Vs is the gradient defined in the 
space E (see Section III). The application of the second law of thermodynamics leads us to 
introduce the positive Onsager matrix C such that 
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^ = A)V E /^ (12) 

This kinetic equation for the ferromagnetic degrees of freedom is actually the simplest 
formulation of the Landau-Lifshitz equation for the magnetization M (see Section II below). 

By analogy with the thermoelectric effect, it is tempting to formally introduce a gradient 
of temperature in the corresponding configuration space. We expect then the existence of a 
supplementary force acting on the magnetization: 



(13) 



J F \ £ S F 
J F I \ rF A / \ V S T 

where S F , A and Tl F are arbitrary Onsager matrices that formally generalises the Seebeck 
coefficient, the thermal conductivity, and the Peltier coefficient in the S space. The question 
is to understand the physical meaning of a temperature gradient in the configuration space 
of the magnetization. This situation is analogous but very different from the case of two 
thermostats of different temperatures localized in two places (in the 3? space). The question 
about the physical signification of a quantity like V^T is not trivial. 

However, the situation described in Eq. ( !T3|) is rather similar to the result obtained at 
the end of the report, providing that the effective temperature gradient Vs^ is replaced by 
the voltage drop due to spin-accumulation: A/i w J ^^dz. A supplementary force is acting 
on the ferromagnet. 

Indeed, let us consider now the system with both spin-dependent electric and ferromag- 
netic dissipation. The ferromagnetic system is not closed, since spins are transfered from 
the electric subsystem to the ferromagnetic subsystem. However, the total system is closed. 
The total internal entropy production allows to access to the kinetic equations of the coupled 
system. The power dissipated is now: 

where J F ot is the total ferromagnetic flux that includes spin transfer and where the vecto- 

— » 6 

rial form of the pumping force A/x is introduced in order to take into account the transverse 
spin-accumulation mechanisms discussed in the litterature related to microscopic theories of 
spin-transfer-torque. 
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Ignoring the electric dissipation (i.e. the two last terms in Eq. (I14p ). the following form 
is obtained for the ferromagnetic system (after some crude simplifications see Section IV): 



7 F \ I £ I 

d tot ' I J~, I 



(15) 



I V E /x J 

where the matrices C, I, I and a are related to measurable experimental parameters. 
In the same way as for the thermoelectric power in Eq. ([6]), the presence of the cross- 
coefficients is justified by the fact that the diffusion of the spin carriers at the interface (e.g. 
s — d relaxation) is contributing to the transport of ferromagnetic moments. 

The consequences of the suplementary term in the expression of the current Jf ot in Eq. (115p 



'tot 

are investigated in terms of a generalized Landau-Lifshitz equation that includes drift and 
diffusion contributions due to spin-transfer [32f (Section IV). The rough arguments presented 
in this introductory section will be developed and detailed in the following sections. 



II. DERIVATION OF LANDAU-LIFSHITZ EQUATION FROM THE CORRE- 
SPONDING POWER DISSIPATION 

In order to treat statistically the time dependence of a unique uniform ferromagnetic 
moment M = M s u r (with radial unit vector u r ) of a fluctuating magnetic nanostructure, 
the ergodic property is used. It allows work with a statistical ensemble of a large number 
of ferromagnetic moments in oriented in the direction {8 ± d9, <p ± dip} of a sphere S of 
radius M s . The density p F (9,ip), defined on the surface of the sphere, is then identified 
with the statistical distribution of ferromagnetic moments. The introduction of the density 
is justified by the nanoscopic size of the magnetic single domain, for which the fluctuations 
play a major role. To that point of view, the system is mesoscopic. Accordingly [47], the 
chemical potential p F takes the general form : 

(j* = kTln{p F ) + V F (16) 

in which the ferromagnetic potential is for instance V F (H, 6) = K sin(9) — M S H cos{6 — 0) 
in the case of a single domain with uniaxial anisotropy of constant K and with an external 
magnetic field H applied at an angle <fi from the anisotropy axis. 



S 



The subscript stands for a closed ferromagnetic system (no source of magnetic mo- 
ments). The corresponding current of magnetic moments Jq is related to the density by the 
conservation law: 

d -§ = -diwJj (17) 

where div-s is the divergence operator defined on the surface of the sphere E. 

The power dissipated by the ferromagnetic system is given by the corresponding internal 

dS F 

entropy production and is given by the product of the generalized flux by the generalized 
force. Assuming a uniform temperature T we have: 

rlQ F 

T M = J ^ VS// ° F (18) 

The application of the second law of thermodynamics dSf /dt > allows the transport 
equation to be deduced by writing the relation that links the generalized flux (the current 
Jq) of the extensive variables under consideration and to the generalized force defined in the 
corresponding space E. Both quantities, flux and forces, are related by the Onsager matrix 
of the transport coefficients C : 

jf = 4W (19) 

The problem is solved as soon as the Onsager matrix is known. In the present case, 
we started from the hypothesis that the magnetic domain was uniform: the modulus of the 
magnetization is conserved. The trajectory of the magnetization (in the configuration space) 
is then confined on the surface of a sphere of radius M s , and the flow is a two component 
vector defined with the unit vectors {u 9 ,ug} of E. Accordingly, the Onsager matrix is a 
2 by 2 matrix defined by four transport coefficients {Lgg, L 9lfi , L^g, L w }. Furthermore, the 
Onsager reciprocity relations impose that L dip = —L v g. 

However, the magnetization is defined by a given axis (unit vector {u r }) in 3D space. 
The choice of the two other vectors is arbitrary, so that L gg = L w . Let us now introduce 
a dimensionless supplementary coefficient a, which is the ratio of the off diagonal to the 
diagonal coefficients: a = Lg^/Lgg. In conclusion, the ferromagnetic kinetic equation is 
defined by two ferromagnetic transport coefficients L 6ip = p$L F and a: 
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£<> --PoL F \ a 1 \ (20) 

-1 a 

On the other hand, the generalized force Vs/z^ , thermodynamically conjugated to the 
magnetization, defines a "generalized" effective magnetic field 

Heff = -V S ^ (21) 



It is a generalization in the sense that this effective field includes the diffusive term [3 



that has first been introduced by Brown in the rotational Fokker-Planck equation 34]. 

The equation Eq. (fT9l) is the well known phenomenological Landau-Lifshitz (LL) equa- 
tion: 



Jo = -PoLf Heff (22) 




Actually, it could be rather surprising to claim that Eq. ( |22|) . that takes the form of 
the Fick's law or thermoelectric laws (with the cross-coefficients), is the "well-known LL 
equation" because the LL equation is the dynamical equation of the ferromagnetic variable 
M = M s u r . However, it is sufficient to rewrite Eq. (122]) in 3D space with re- introducing 
the radial unit vector u r = (1,0,0) of the reference frame {u r , ug, u^}, and recalling that 
the current is the density multiplied by the velocity Jq = p^dur/dt, in order to recover the 
traditional LL equation from ( !22l : 



du r 
~dt 



jw r x H e fj — au r x (u r . x H e ff^ (23) 



Furthermore, it is well-known that LL equation is equivalent to the following Gilbert 35 1 
equation, that allows the damping coefficient t] to be defined: 

-^xr^-^M.-j (24) 

where V is the gyromagnetic ratio. The equivalence between the two equations defines 
the coefficients a and Lp has a function of the coefficients rj and T. a is the dimentionless 
damping coefficent: 
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a = rjTM s 

and Lp is defined by the relation 



(25) 



M s (1 + a 2 



(26) 



The corresponding Fokker-Planck stochastic equation first derived by Brown 34| is ob- 
tained directly by inserting Eq. ( 122"]) into Eq.((T] 



dt 



-V e £ Ve//o 



(27) 



Using Eq. ( |T6|) . Eq. ( |25|) . Eq. ( |26l) . and the explicit expression of the Laplacian V| in 
spherical coordinates, the Fokker-Planck equation reads: 



dpp _ L F d 
dt ~ sin 6 06 



+ 



sin 6 



d 



dV F 1 dV 



06 
dV F 



sin 6 d<j) 
a dV \ 



a dpf 



06 



(28) 



sin 6 d<p \ \ 86 ' sm.6 d<p ) ru sin^ 
The driving force responsible for the magnetization dynamics is distributed between drift 
and diffusion terms for the probability distribution. The equilibrium solution of the equation 
is the Boltzmann distribution, as it was assumed in expression Eq. (1161) for the definition 
the chemical potential . 

Since the equation depends on the determinist potential V F (that contains the energy 
due to the external magnetic field, the magnetocrystaline anisotropy, dipolar energy, etc), 
it is non linear. Only a few simple configurations can find an analytical solution for the 

rn, 

non-equilibrium statistical distribution p F (6,<p,t) |36|. This is typically the case for linear 
expansions near equilibrium states in the context of ferromagnetic resonnance, or for the 
Neel-Brown activation process 34|, |36( at long time scales. Eq. (T25|) will be extended to 



spin-transfer contributions introduced in Section III, after the study of spin-accumulation 
effects below. 
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III. DERIVATION OF SPIN- ACCUMULATION FROM THE CORRESPONDING 
POWER DISSIPATION 



In this section, we focus on the electric transport only (we forget the role of the fer- 
romagnetic variable and the existence of the S space). The corresponding electric wire is 
defined along the z axis, with a section unity: the relevant configuration space is the one 
dimentional real space 9ft. However, in order to take into account the spin-dependent electric 
current, the two channel model is introduced. Beyond the diffusive two chanel model ap- 
proach, transverse spin accumulation parameters are also introduced, in order to take into 
account the spin precession in quasi-ballistic regime near the interface. 



A. two channel relaxation and spin-diffusion 

The system is composed of two populations of conduction electrons with a relaxation 
process that allows the electrons of one population relax into the other. The difference 
between the two populations is introduced through an internal degree of freedom, and the 
relaxation process occurres within the space defined by this internal degree of freedom. In- 
side the bulk, the relaxation from one channel to the other is compensated by the opposit 
relaxation: the electronic populations are maintained at equilibrium. However, the presence 
of an interface with inhomogeneous transport parameters puts out-of-equilibrium the elec- 
tronic populations. At steady state, a diffusion process (of the spin-density) in the 9ft space 
occurres that compensates the forced relaxation defined in the internal space. This diffusion 
process is called spin-accumulation in the case of spin-dependent transport. Formally, the 
two-channel model consisted of defining a bi-valuated internal variable for the transport 
parameters, that takes the values a or 7. Typically, the values of the internal variables are 
a =t and 7 =\, for spin-flip scattering, or a = s (f l) an d 7 = d (4-) for spin-dependent 



s — d scattering 16, 26] ( the three channel model in which the internal variable takes the 

n 

values {s "|, s |, d \) is developed in reference [16J). Accordingly, the local electrochemical 
potentials are defined by [f a and /i^, and the electric currents generated in each channel is 
noted {Ja,Jy} 

The conservation laws write: 
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where n a and n 7 are the densities of charge carriers in the channels {0,7}, and the spin- 
dependent relaxation is taken into account by the flux St. This is the velocity of the reaction 
(or relaxation of the spin-dependent internal variable) that transforms a conduction electron 
a into the conduction electron 7. This generalized flux defines a "spin current" (density times 
velocity) in the configuration space of the internal variable (somehow related to S: see next 
section). Note however that in the litterature the term "spin current" is devoted to the 
spin-polarized electric current 5J e = — J 7 defined in the real space 3?. 

The power dissipated by the electric system is given by the corresponding internal entropy 
variation, i.e. by the product of the currents by the electric fields: 

where we introduced the difference of the chemical potentials A/i e = fi% — [13|, Il4j . 
the application of the second law of thermodynamics leads to the kinetic equations, after 
introducing the transport coefficients: the conductivities a a , cx 7 , and the Onsager coefficient 
L, such that: 

je erg dfMg 

J a e dz 

J l ~ e dz K J 

* = LAfi e 

where the two first equations are Ohm's law applied to each channels. The effect of the 
electric charge relaxation is described in reference [16]. The Onsager coefficient L is shown to 
be inversely proportional to the electronic relaxation times t q< _ >7 . The total electric current 
is spin-independent: 

Jo = J £ a + Jy = ~\y z + 0-7/4) ( 32 ) 

However, it is not possible to measure separately the different conduction channels, since 
any realistic electric contact short cuts the two channels. What is measured is necessarily 
the usual Ohm's law, Jq = — 2<To^, that imposes the reference electric potential ( to be 
introduced, together with the mean conductivity a = (cr Q , + a 7 )/2. The potential ( is hence: 
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eC 



Co 



(33) 



The reference configuration is defined by the two channels collapsing to a unique con- 
duction channel (e.g. parallel magnetization of a junction of two identical ferromagnetic 
layers: Ajjf^ (0) = ). The non-equilibrium (A/i e (0) 7^ 0) contribution of the junction to 
the resistance, R ne , is calculated through the relation: 



so that 



JleK 



d_ 

dz 



(/4 ~ e((z))dz 



d_ 

dz 



(/£ - e({z))dz 



(34) 



B 



a a — a 7 dA^ 



R™ = -ir ""—^7 dz (35) 

Jo e J A ^0 dz 

where the measurement points A and B are located far enough in each side of the interface 
(inside the bulk) so that Afi e (A) = Afi e (B) = 0. The integral in Eqs. fl34|) is performed over 
the regular part of the function only (across the interface ( and <Tj are discontinuous at this 
scale): this resistance is proportional to the discontinuity at the interface. It is convenient 
to describe the conductivity asymmetry by the parameter (3 such that a a = a (l + (3) and 
cr 7 = er (l — (3). On the other hand, the spin-polarized electric current is 5J e = — Jf;. 
With these new variables, Eq.f l35|) rewrites: 



R ri 



and the Onsager matrix reads: 



B 



2f3_ 



dAf/ 

dz 



dz 



(36) 



J' 



\ / -1 dAp* \ 
e dz 
-1 dul 



(TO /?cr 

^0 a ^ d -§ (37) 

L j \ Afi j 
The system of equations Eq. (1371) allows the diffusion equation for A/x(z) to be derived 
for the stationary conditions df^/dz = and 85 ~J /dz = —2^: 



V 



d 2 A/i Afi 

dz 2 



L diff 



(38) 



where l di ff 
configurations 



~0 



16 



t£z m • The resistance R ne can then be calculated for each specific device 
26). 
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The two channel approximation with the internal degree of freedom that takes the value 
{t, 1} describes the consequences of the spin-flip scattering. However, the model is also 



sufficient for the description of spin-dependent s — d relaxation [16l . l26j . where the d band 
is full for the majority spins t- Indeed, the s — d relaxation with spin-flip (from s t to d \.) 
has a very small probability to occurre, compared to the relaxation without spin-flip (from 
s l to d 1). As a consequence, the two channel model can also be used, and it also leads to 
a redistribution of the spin populations at the interfaces, i.e. to spin-accumulation. In both 
cases, the spin-accumulation is described by the function Afi(z), which is solution of is the 
the diffusive equation Eq. ( I38lh 



B. Quasi-ballistic effect and transverse spin-accumulation 

However, the description proposed above with a spin-dependent internal variable that 
takes the two values {a, 7} is not able to take into account the precession of the spins 
occuring in a magnetic field. If the precession contribution is not relevant in the case of the 
processes that lead to giant magnetoresistance (because the mean values are averaged out 
over the spin- diffusion length) this is no longer the case in a quasi ballistic regime near the 
interface. 

In order to take into account quasi-ballistic effects near the interface (i.e. sub-nanometric 
scales in metalic devices), the two-channel model has been recently generalized to trans- 



verse spin-accumulation in the context of spin-transfer-torque investigations |37H42|. The 
transverse spin-accumulation is introduced with the corresponding current 5J]_ and the 
corresponding chemical potential Afi e ± . Transverse means here that the spin density is con- 
sidered in the plan perpendicular to the quantification axis £ that defines the spin up and 
spin down in the two channel-model. 

The coefficient a± can also be defined through the corresponding diffusion coefficient 



D± = where n± is the density of transverse spins |43j . It is then also possible to define 
a "pseudo" spin-diffusion process in the case of spin-decoherence. Note however that the two 
potentials A/i e and Afi e ± are defined at very different length scales and it is necessary to refer 
to quantum approaches in order to understand the physical signification of the transverse 



parameters 



44 



45] . The corresponding contribution to the power dissipated is 
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T 



dS e t 



-6J1 



dAfi e ± 



dt " M edz 

Puting all together we have the following Onsager relations for the electric system: 



(39) 



( 5.P\ 
Jo 

V * j 



f (T f3a \ 

f3a <t 

<r ± 

\ L J 



/ -1 dAfi e \ 

e dz 

-1 d/^g 
e cfe 

-l d a m i 



e 3,2 



(40) 



IV. DERIVATION OF SPIN-TRANSFER DUE TO LONGITUDINAL AND 
TRANSVERSE SPIN- ACCUMULATION 



In usual experimental configurations for spin-transfer, an electric current is injected in 
a circuit that includes a ferromagnet (in series or in non-local configuration 48[) and the 
magnetoresistance, i.e. the potential drop (of the form — § £f ^^dz Eq. (|36|) ) allows 
the magnetization states to be measured. The effect of strong electric currents on the 
magnetization states can then be observed. In such a configuration, the two sub-systems 
described in the sections above exchange magnetic moments at the junctions and both are 
open systems. 

In order to describe the dynamics of the ferromagnetic degrees of freedom (following step 
by step the method presented in Section II), we have to deal with a closed system. The 
system of interest is now the ferromagnetic system that includes spin-accumulation effects 
at the junctions. This total ferromagnetic system is such that the density of ferromagnetic 
moments pf ot and the total ferromagnetic flux J[ ot are related by the conservation law: 
dpfot/dt = -div s Jf ot . 

The initial configuration space of magnetic moments is then extended to ID real space 
parametrized by the internal variable E®K a7 . The important point here is that the internal 
variable is spin-dependent, and related to the ferromagnetic space S (e.g. through s — d 
relaxation and the corresponding spin-accumulation). This accounts for the coupling, i.e. 
the transfer, of magnetic moment between the two sub-systems. 

The dissipation is given by the internal power dissipated in the total system TdSi/dt : 
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dSi -?F <n ,,F sr re 9A^J, e <9A/1^ -, e d\l% ,f rA „e ( A -x\ 

T^=j to ,V s/ , -tfJ._- ( yj ± ._-j _ + *A A4 (41) 

Where the first term in the right hand side is the power dissipated by the total ferro- 
magnetic sub-system (including the ferromagnetic contribution due to spin-transfer), the 
two following terms are the power dissipated by spin-dependent electric transport, and the 
fourth term is the spin-independent Joule heating. The last term is the power dissipated by 
spin-flip or s — d relaxation. 

In Eq. f|4Tj) . the vectors are defined on the sphere £ with the help of two angles 9 and 
ip. The total ferromagnetic current jf ot = j[ Q fue + jtot^ includes the contribution due to 
spin-accumulation mechanisms. The chemical potential pF accounts for the energy of a 
ferromagnetic layer. On the other hand, the system is contacted to electric reservoirs with 
the electric currents and the corresponding chemical potentials. Applying the second law of 
thermodynamics, we obtain the general Onsager relations: 
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(42) 



All coefficients were defined in the previous sections, except the new cross-coefficients 
{lip V , l<pe, he, hip), introduced in this model as spin-transfer coefficients, related to the experi- 
mental parameters. The coefficients {li} i= ^ e ^ are deduced from the coefficients {/j} through 
the Onsager reciprocity relations. 

The total ferromagnetic current can be written after integrating over the volume v of 
the ferromagnetic layer of section unity and the spin accumulation zone. This volume is 
such that v = dz, where z = A and z = B are two sections close to the interface but 
far enough with respect to the diffusion length Idiff- We assume here that hiff is much 
smaller than the width of the ferromagnetic layer in order to simplify the calculation: the 
volume of the ferromagnetet is identified as v. Let us define X as the correction due to the 
spin-transfer deduced from the two first equations of the matrix equation Eq. (j4~2]) . after 
integrating over the volume v : 
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vJZ t = vC Vrf + X (43) 
where Co is the matrix defined in Eq. ( |20l ). 

The assumption of constant modulus of the magnetization imposes that X is confined on 
the surface of the sphere S. The Helmoltz decomposition theorem can then be applied: the 
vector X can be decomposed in a unique way with the introduction of the two potentials x 



and $ (i.e. a potential vector) such that 49]: 



X = u r x V E $ + VsX (44) 

where the first term is divergenceless and the second term is curless (i. e. non conserva- 
tive). The method used here is hence not equivalent to that of adding a spin-transfer source 
term in the conservation equation of the time variation of pf ot . The two potentials will be 
described in more details below. 

The total correction to the Landau-Lifshitz-Gilbert equation writes: 

v Jf ot = w£ V s /it + i x V E $ (45) 

where the electro-spin chemical potential that describes the ferromagnetic system with 
the addition of the spin-transfer contribution writes 

Ht = kTln(p F ) + V F + ^- (46) 

However, in Eq. (T461) . the density p F is no longer relevant because the ferromagnetic 
system alone is an open system, and only the total density p F ot is defined. The canonical 
form of the chemical potential of the total system that contains the total density p to t is: 
/if = kTln(pf ot ) + V F . The total density is deduced with identifying with Eq. ( )46l) p\ 



F 
of 



The generalized LLG takes the form: 



d,U r 



7 



dt MJl + a 2 ) 



u r x I H e f f + 



vp F L F 







j — au r x 


U r X | 



l eff 



vp F Lpa 



(47) 



This is the main result of this work. 
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The corresponding Fokker-Planck equation is obtained by inserting the expression of J to t 
into the conservation equation: dp[ ot /dt = —div-sJ^. The study of the resulting stochastic 
equation is however beyond the aim of this paper. 

In order to give an expression of the two potential-energy terms {x, 3>}, we will make the 



following assumption. According to previous discussions [26|, |32|, |46j based on the separation 
of the typical relaxation time scales involved during the ferromagnetic processes, the usual 
spin-accumulation due to spin-dependent relaxation Afi e is coupled to the relaxation of the 
magnetization jf£ only (because Idiff defines a mesoscopic variable that scales with the 
magnetization). On the other hand, we assume that the transverse spin accumulation is 
coupled to the precession only (i.e. acting at subnanosecond time-scale). 

In this case and after integrating over the volume v the ferromagnetic current writes : 

„ TFd — 7; jFO _i_ d\ l_d® 

UJ tot — UJ do sine dip (4gN 

j F f _ j F f I a* 

VJ tot — VJ ' 90 

On one hand the potential x is directly associated to the voltage drop due to spin- 
dependent relaxation (Eq. 136]) : 



dz= _mZL B UMK (49) 

09 e J A dz 2(3 K J 

and this expression can be generalized to specific device configurations. The potential 
energy x, function of the magnetic coordinates, can be measured in the context of two-level- 



fluctuation experiments performed on individual magnetic nanostructures [32|, |46 |. 

On the other hand, the potential $ is associated to the discontinuity of the transverse 
spin accumulation and is responsible for the spin-transfer torque: 

lysine rB 9A/U r 

e J A dz f^fl"! 

e J A dz 

In these expressions, and in analogy with spin-accumulation due to the spin-dependent 
relaxation, the voltage drop ^^dz = J^R^f 18 is also able to define a non-equilibrium 
interface magnetoresistance R^J 13 for transverse spin-accumulation. Experimentally, the 
potential $ is identified to the so called Slonczewski term used in the context of resonance 
experiments (FMR) performed in the GHz range. 
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V. CONCLUSION 



In order to describe spin-transfer effect, spin-accumulation has been taken into account 
explicitly in the dynamical equation of the macroscopic ferromagnetic degrees of freedom. 
This dissipative coupling was described in terms of Onsager cross-coefficients U appearing 
in the Onsager matrix in Eq. (1421) . 

All the other terms appearing in the equation Eq. fj42|) (the other transport coefficients 
of the matrix, the generalized flux, and the conjugated generalized forces) have first been 
defined independently in the two preceding sections. In the case of the well-known dynamics 
of the ferromagnetic order parameter (supposed uniform), the approach proposed allows us to 
define the transport coefficients, the ferromagnetic current and the ferromagnetic generalized 
force from the expression of the entropy production and the conservation equations. The 
dynamics of the magnetization is summarized in the Onsager equation Eq. (T22|) . which 
is the simplest form of the well-known Landau-Lifshitz equation. On the other, the spin- 
dependent electronic relaxations (spin-dependent s—d relaxation or spin-flip relaxation) were 
treated on an equal footing in the context of the two channel model of electric conductivity. 
The resulting kinetic equations are also summarized by a Onsager equation Eq (l57|) with 
the relevant flux and forces, which are also a simple form of well known kinetic equations 
(e.g. that derived from the Valet-Fert model). The corresponding transport coefficients 
and forces (/3 , cr , SJ e , A/i) can be measured through the giant magnetoresistance and 
related effects. Due to quasi-ballistic precession of the spins of the conduction electrons, 
it is furthermore necessary to generalize the spin-accumulation effect to "transverse spin- 
accumulation" (according to recent reports on spin-transfer torque). 

Due to the two forms of spin-accumulation mechanisms, there are also two forms of cou- 
pling, namelly spin-accumulation coupling (due to spin-dependent relaxations) and trans- 
verse spin-transfer-torque coupling due to quasi-ballistic spin precession. The model pro- 
poses a method able to formalize this coupling, and to deduce the consequences in terms of 
Landau-Lifshitz equation. In both cases, the coupling between the spin of the conduction 
electrons and the ferromagnetic parameters are introduce through the four phenomenolog- 
ical Onsager cross-coefficients Uj. The generalization of the Landau-Lifshitz equation to 
these contributions is performed with two measurable potentials (functions of the 

magnetic coordinates) defined in Eqs. ( l45l) and ( 1461 ). The potential x is associated to 
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the spin-accumulation generated by the spin-dependent relaxation, and the potential <E> is 
associated to the conservation of the transverse moments (spin-transfer-torque). The two 
functions are experimentally accessible. 
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